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Abstract 

The cross section for e^e~ 7r~^7r~ J/i/j between 3.8 and 5.5 GeV/c^ is measured using a 548 fb~^ 
data sample collected on or near the T(4S') resonance with the Belle detector at KEKB. A peak 
near 4.25 GeV/c^, corresponding to the so called y(4260), is observed. In addition, there is another 
cluster of events at around 4.05 GeV/c^. A fit using two interfering Breit-Wigner shapes describes 
the data better than one that uses only the y(4260), especially for the lower mass side of the 
4.25 GeV enhancement. 

PACS numbers: 14.40. Gx, 13.25.Gv, 13.66.Bc 
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In a recent study of initial state radiation (ISR) events of the type, e~^e~ — >■ 
T/SiiTr+vr" J/-?/;, the BaBar Collaboration observed an accumulation of events near 
4.26 GeV/c^ in the tt^tt" J invariant mass distribution and attributed it to a possible new 
resonance that they dubbed the y(4260) This observation was confirmed by the CLEO 
experiment using a similar technique with a data sample collected at the T(45') peak 0]. 
The CLEO Collaboration also collected a 13.2 pb~^ data sample at ^/s = 4.26 GeV, and 
reported signals for 7r~^7r~ J/ip, tt^ti^ J /ip, and K^K~ J/ip with cross sections that are signif- 
icantly higher than those measured at other nearby energies [3|. 

Since the y(4260) resonance is produced via e^e~ annihilation accompanied with initial 
state radiation, its J^'^^ = 1 . However, the properties of the observed peak are rather 
different from those of other known J^*^ = 1 charmonium states in the same mass range, 
such as ?/'(4040), ip^ilGO), and ipi^Mld). Since it is well above the DD threshold, it is 
expected to decay predominantly into final states. The partial width for the tt-kJ /ip 

final state is expected to be a small fraction of the total. In fact, the 1^(4260) shows an 
unusually strong coupling to the TnTj/ip final state while no significant enhancement is 
observed in D^*^D^*^ final states In a fit to the total hadronic cross sections measured 
by the BES experiment [s], 0] for ^/s between 3.7 and 5.0 GeV, Mo et al. set an upper 
limit on T^+e- for the F(4260) to be less than 580 eV at 90% confidence level (C.L.) 0. 
This implies that its branching fraction to -KirJ/ip is greater than 1.3% at 90% C.L. These 
properties have triggered many models to explain the y (4260) as an exotic state, such as a 
four-quark state, a molecular state, or a quark-gluon hybrid |8|. 

In the analysis reported here, we use a 548 fb~^ data sample collected with the Belle 
detector operating at the KEKB asymmetric-energy e~^e~ (3.5 on 8 GeV) collider 10] to 



investigate the 7r~^7r~ J/ ijj final state produced via ISR. About 90% of the data were collected 
at the T(4S') resonance {^/s = 10.58 GeV), and about 10% were taken at a center-of-mass 
(CM) energy that is 60 MeV below the T(4S') peak. The measurement in this Letter uses 
an irnproved efficiency for detecting ISR events, and supersedes the preliminary results in 
Ref. [ll|, which confirmed the structure near 4.26 GeV/c^. 

For Monte Carlo (MC) simulations of the ISR process, we generate signal events with 
the PHOKHARA program [13]. In this program, after one or two photons are emitted, the 
lower energy e^e~ pair forms a resonance X that subsequently decays to n^n^J/ip with 
the J/ip decaying either to e~^e~ or fi'^fi~. In the X ti'^tx" J/ip generation, we use pure 
S*- waves between the tttt system and the J/'0, as well as between the ir^ and vr"; this is in 
agreement with the experimental results [l|, [l3| . The -k^ti' invariant mass distributions are 
generated according to phase space. For ip{2S) — > 7r~^7r~ J/ip, which we use as a calibration 



process, we use the decay properties that have been measured with high precision [13[ . 

For candidate events, we require the number of charged tracks to be four and net charge 
to be zero. For these tracks, the impact parameters perpendicular to and along the beam 
direction with respect to the interaction point are required to be less than 0.5 and 4 cm, 
respectively, and transverse momentum is restricted to be higher than 0.1 GeV/c. For 
each charged track, information from different detector subsystems is combined to form a 
likelihood for each particle species (i), Li IJ]. Tracks with 11,^ = Ck+c ^ identified 



as pions with an efficiency of about 95% for the tracks of interest. Similar likelihood ratios 
are formed for electron and muon identification. For electrons from J/ip ^ e+e~, one 
track should have TZ^ > 0.95 and the other TZ^ > 0.05; for muons from J/ip ^ 
at least one track is required to have TZ^ > 0.95; in cases where one of the tracks has 
no muon identification (ID) information, the polar angles of the two muon tracks in the 
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TT+vr"/!"'"/!" CM system are required to satisfy | cos6'^| < 0.7 based on a comparison between 
data and MC simulation. Lepton ID efficiency is about 90% for J/ip ^ e^e~ and 87% for 
J/ip ^ H~^IJ'~- Events with 7-conversions are removed by requiring TZ^ < 0.75 for the tt+tt" 
tracks. For the J/ip ^ e^e~ mode, 7-conversion events are further removed by requiring 
the vr+vr" invariant mass to be greater than 0.35 GeV/c^. 

The detection of the ISR photon is not required, instead, we identify ISR events by the 
requirement \M^^^\ < 2.0 (GeV/c^)^, where M^^^ is the square of the mass that is recoihng 
against the four charged tracks. 

Clear J/ip signals are observed in both decay modes. We define a. J ftp signal region as 
3.06 GeV/c^ < m^+^- < 3.14 GeV/c^ (the mass resolution is about 17 MeV/c^), and J/^jj 
mass sidebands as mi+i- G [2.91,3.03] GeV/c^ or m^+^- G [3.17,3.29] GeV/c^; the latter are 
three times as wide as the signal region. 

Figure [T] shows the invariant mass [l5| distribution after the above selection, 

together with the background estimated from the J / ip mass sidebands. In addition to a huge 
ip{2S) signal, there is a clear enhancement at 4.25 GeV/c^ similar to that observed by the 
BaBar Collaboration [l[. In addition, there is a clustering of events around 4.05 GeV/c^ that 
is significantly above the background level. It is evident in the figure that the background 
estimated from the J/ip sidebands agrees well with the level of the selected events in the 
high 7r+7r~£+£~ invariant mass region. A study of events in the |M^^g^| > 1 (GeV/c^)^ region, 
which is depleted in signal events, supports this conclusion. The backgrounds not in the 
sidebands, including: (1) vr+vr" J/-?/;, with J/ip decays into final states other than lepton 
pairs; (2) XJ/ip, with X not being vr+vr", such as K^K~ and vr+vr"-/?^, are found from MC 
simulation to be less than one event per 20 MeV/c^ bin at 90% C.L. according to the CLEO 
measurements jsf and are neglected. The production of ti^tt' J /ip from non-ISR processes, 
such as e~^e~ — > 77*7* —>■ 'jp^J/i/j, is computed to be small [l6| and is neglected. 
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FIG. 1: Invariant mass distribution of 7T~^TT~i~^i~ . The blank histograms represent the selected 
data and the shaded histograms are the normalized sidebands. The inset shows the distribution 
with a logarithmic vertical scale. 
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The data points in Figs. M^a) and (b) show the background-subtracted M^^^ distribution 
and the polar angle distribution of the tt^h' J /ip system in the e^e~ CM system for the 
selected tt^'k'J/iI) events with invariant mass between 3.8 and 4.6 GeV/c^. The data agree 
well with the MC simulation, indicating that the signal events are produced via I SR. 
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FIG. 2: M^g^ distribution (a) and the polar angle distribution of the tt^tt^J/i/j system in the e'^e~ 
CM frame (b) for the selected Tr~^Tr~J/ip events with invariant masses between 3.8 and 4.6 GeV/c^. 
The background from J/tp mass sidebands has been subtracted, and the selection criterion applied 
to the M^^j, has been relaxed in (a). The points with error bars are data, compared with MC 
simulation (solid histograms). 




We estimate the signal significance of the clusters at 4.05 GeV/c^ and 4.25 GeV/c^ by 
comparing the numbers of signal events (number of observed events in the J/ signal window 
minus the number of J/^/'-sideband-estimated background events) with their statistical un- 
certainties. For events with m^+^-e+e- e [3.80,4.15] GeV/c^, we have n''s(4.05) = 120 ± 14, 
which is more than 80" from zero assuming a Gaussian error; while for events with 
m^+^-i+i- e [4.15,4.60] GeV/c^, we have n^'g(4.25) = 324 ± 21, which is more than 15a 
from zero. 

The e^e~ tt^tt^J/iP cross section for each TX^-K^J/ip mass bin is computed with 

(Ji = i;-c-B(j/i,^e+£-) ' where n° , , Ei, and Ci are the number of events observed in data, 
the number of background events determined from the J/ip sidebands, the efficiency, and 
the effective luminosity |l7| in the 2-th Tr~^7r~J/ip mass bin, respectively; B{J/ip = 
11.87% is taken from Ref. [l8|. The resulting cross sections are shown in Fig. [3l where the 
error bars indicate the combined statistical errors of the signal plus background events. Our 
measurement at 4.26 GeV/c^ agrees well with BaBar's and CLEO's results [l|, llj. 

The sources of the systematic errors for the cross section measurement are listed in Ta- 
ble [H The particle ID uncertainty, measured using the ip{2S) events in the same data sample, 
is 3.0%; the uncertainty in the tracking efficiency for tracks with angles and momenta charac- 
teristic of signal events is about 1%/track, and is additive; efficiency uncertainties associated 
with the J / ip mass and M^^^ requirements are also determined from a study of the very pure 
ip{2S) —>■ TT^TT~ J/ip event sample. In this study we find that the detection efficiency is lower 
than that inferred from the MC simulation by (2.5 ±0.4)%. A correction factor is applied to 
the final results and 0.4% is included in the systematic error. Belle measures the luminosity 
with a precision of 1.4% using wide angle Bhabha events, and the uncertainty of the ISR 
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FIG. 3: The measured e"*"e tt^tt J ftp cross section for CM energies between 3.8 and 5.5 GeV. 
The errors are statistical only. 



photon radiator is 0.1% [17 



12| generator is 



The main uncertainty of the PHOKHARA 
due to the modelhng of the tt+tt" mass spectrum. Figure H] shows the vr+vr" invariant mass 
distributions of events for three m^+^-j/^ regions, [3.8,4.2], [4.2,4.4], and [4.4,4.6] (unit 
in GeV/c^). The tt+tt" invariant mass distribution for events around 4.25 GeV/c^ differs 
significantly from phase space; for other energy ranges the agreement with phase space is 
better. Simulations with modified vt+tt" invariant mass distributions yield efficiencies that 
are higher by 2-5% for nij^+j^-j/^ below 4.4 GeV/c^. This is not corrected for in the analysis, 
but is taken as the systematic error (conservatively assigned as 5%) for all 7r~^7r~J/ip mass 
values. The selected events have four charged tracks and 16-25% of them have a detected 
high energy ISR photon. According to the MC simulation, the trigger efficiency for these 
events is around 98%, with an uncertainty that is smaller than 1%. The uncertainty of 
B{J/il) — * = B{J/tp —> e+e~) + B{J/ip is taken as 1% by hnearly adding the 

errors of the world averages for the e+e~ and fi~^fi~ modes [l^. Finally the MC statistical 
error on the efficiency is 1%. We assume all the sources are independent and add them in 
quadrature, resulting in a total systematic error on the cross section of 7.5%. 

As a validation of our analysis, we measure the ip{2S) cross section with the same selection 
criteria. Here 15,444 ip{2S) events survive the selection and the MC-determined detection 
efficiency is 5.13%. This corresponds to ct(?/'(2^)) = (15.42 ± 0.12 ± 0.89) pb at the T(4^) 
resonance or r{ilj{2S) — e~^e~) = (2.54 ±0.02 ±0.15) keV, where the first error is statistical 
and the second systematic. This measurement agrees well with the world average value of 
(2.48 ± 0.06) keV [l8|. The ipi^S) mass determined from the data indicates the 7r~^7r~i^i~ 
invariant mass is measured with a precision of ±0.6 MeV/c^. 

An unbinned maximum likelihood fit is applied to the 7r'^7T~i~^i~ mass spectrum in Fig. [H 
Here the theoretical shape is multiplied by the efficiency and effective luminosity, which are 
functions of the 7T^7r~i^i~ invariant mass. Since there are two clusters of events in the mass 
distribution, we fit it with two coherent Breit-Wigner (BW) resonance functions (-R1, R2) 
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TABLE I: Systematic errors in the cross section measurement. They are common for all data 
points. 



LAX V-'VJ 


Relative error iTni 


Particle ID 


3.0 


Tracking 


4 


J /if) mass and M^g^, selection 


0.4 


Integrated luminosity 


1.4 


m^+^- distribution 


5 


Trigger efficiency 


1 


Branching fractions 


1 


MC statistics 


1 


Sum in quadrature 


7.5 




FIG. 4: The 7r"'"7r invariant mass distribution of events for different tt^tt J/ip mass regions, 
(a): m^+^-j/^ G [3.8,4.2] GeV/c^, (b): m^+^-j/^, £ [4.2,4.4] GeV/c^, and (c): m^+^-j/^ G 
[4.4,4.6] GeV/c^. The points with errors bars are pure signal events, the histograms are MC 
simulations made using phase space distributions. 

assuming there is no continuum production of e'^e^ TT^n^J/ip. In the fit, the background 
term is fixed at the level obtained from a linear fit to the sideband data, contributions 
from the ip{2S) and '?/;(3770) resonance tails (added incoherently) are estimated using world 
average values for their parameters [l^ and fixed, the widths of the resonances are assumed 
to be constant. A three-body decay phase space factor is applied. The MC-determined mass 
resolution is less than 5 MeV/c^ over the full mass range. This is small compared to the 
widths of the resonances in our study and is ignored. 

Figure [5] shows the fit results; there are two solutions with equally good fit quality. The 
masses and widths of the resonances are the same for both solutions; the partial widths to 
e~^e~ and the relative phase between them are different (see Table [TTl) [3| • The interference 
is constructive for one solution and destructive for the other. The systematic errors come 
from the absolute mass scale, the detection efficiency, the background estimation, the phase 
space factor, and the parametrization of the resonances. The quality of the fit assessed from 
the binned distribution of Fig. [5l is /''^df = 81/78, corresponding to a C.L. of 38%. The 
statistical significance of the structure around 4.05 GeV/c^ is estimated to be 7.4cr from the 
change in likelihood value when the BW representing it is removed from the fit. Although 
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the mass of the first resonance is close to that of the iplAOAO), the fitted width is much 
wider than its world average [l8| value (80± 10 MeV/c^). The mass of the second resonance 
is higher than that of the ip{4:160). Changes of resonance parameters that occur when we 
fit with a coherent ip{2S) tail, a coherent or incoherent non-resonance term, an energy- 
dependent total width, or a cascade two-body phase-space factor, dominate the systematic 
errors listed in Table Hit the significance of the Rl signal is greater than 5cr in all of the fitting 
scenarios that are considered. If we use the same functional form as BaBar (a single BW with 
an incoherent second-order polynomial background term) we find M = 4263 ± 6 MeV/c^, 
Ftot = 126 ± 18 MeV/c^, and Bin+n-J/tp) ■ T^+e- = 9.7 ± 1.1 eV/c^, consistent with their 
results [l|. 




FIG. 5: Fit to the tt^ tt" J / ip mass spectrum with two coherent resonances. The curves show the 
best fit and the contribution from each component. The dashed curves are for solution I, and the 
dot-dashed curves for solution II. The histogram shows the scaled sideband distribution. 



In summary, the e~^e~ tt^tt^J/iP cross section is measured for the CM energy range 
^/s = 3.8 GeV to 5.5 GeV. There are two significant enhancements: one near 4.25 GeV, 
consistent with the results of Refs. [H and and another near 4.05 GeV, which has not 
previously been observed. We note that these enhancements are close to D thresh- 
olds, where coupled-channel effects and rescattering may affect the cross section If 
we nevertheless represent the cross section using interfering BW terms, a second term (in 
addition to the y(4260)) substantially improves the fit. In particular, the lower-mass side 
of the 4.25 GeV enhancement is better reproduced. The parameters that are obtained from 
this two-term fit do not correspond to those of any of the excited ip states currently listed 



in Refs. 18 and 21 
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TABLE II: Fit results of the vr+vr J /ip invariant mass spectrum. The first errors are statistical and 
the second systematic. M, Ftoti and B ■ Fg+g- are the mass (in MeV/c^), total width (in MeV/c^), 
product of the branching fraction to 7r~^7r~J/ijj and the e~^e~ partial width (in eV/c^), respectively. 
(j) is the relative phase between the two resonances (in degrees). 



Parameters Solution I Solution II 





M{Rl) 


4008 ± 40+^^^ 




rtot(i?i) 


226 ± 44 ± 87 


B 


• Fg+g-(m) 


5.0±1.4t^;^ 12.4 ± 2.4l}^;L^ 




M{R2) 


4247 ± 12l^^ 




rtot(i?2) 


108 ± 19 ± 10 


B 


•Fg+g-(i?2) 




6.0±1.2t^:^ 20.6±2.3l?;^ 
12±29t^g -lll±7j:^^ 
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